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On the coherent effect of vacuum fluctuations on driven atoms 


M. Donaire and A. Lambreclit 

Laboratoire Kastler Brossel, UPMC-Sorbonnes Universites, ENS-PSL Research Universities, 
CNRS , College de France, place Jussieu, F-75252 Paris, France 
(Dated: 10 April 2015) 

We study the coherent effect of the Casimir-Polder interaction on the oscillations of two-photon 
driven atoms. We find that, for oscillations between two degenerate states in lambda-configuration, 
shifts on the Rabi frequency may be induced by non-additive level shifts. For oscillations between 
two Rydberg states in ladder-configuration, shifts on the Rabi frequency may be induced by the 
effective renormalization of the laser interaction. 
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I. INTRODUCTION 

The interaction of a neutral atom with a material sur¬ 
face is a problem profusely addressed in the literature 
In most of the approaches the atom is taken in 
a stationary state w.r.t. the time of observation. At 
zero temperature and in the electric dipole approxima¬ 
tion, the atom undergoes a series of virtual El transitions 
to upper levels. It is the coupling of the charges of the 
atom and the currents on the surface to the quantum EM 
field that induces the correlation between their transtient 
dipole moments, giving rise to a non-vanishing interac¬ 
tion. The lifetime of the virtual atomic transitions is 
very short in comparison to ordinary observation times 
and thus, the use of stationary quantum perturbation 
theory is well justified for the calculation of this interac¬ 
tion [2j. For distances greater than the relevant atomic 
transition wavelengths this interaction is referred to as 
retarded Casimir-Polder (CP) interaction , while for dis¬ 
tances much shorter than those wavelengths it is referred 
to as van der Waals or non-retarded CP interaction. 

When atoms close to dielectric surfaces are driven un¬ 
der the action of external sources, virtual and actual tran¬ 
sitions mix up with each other in the time evolution of 
the atomic wave function. From a practical perspective, 
the effect of the CP interaction on the dynamics of driven 
atoms is of great importance in hybrid quantum systems 
-eg. Ref.jgJ. At first sight, under quasi-stationary condi¬ 
tions, that effect reduces to an additive level shift on the 
atomic eigenstates, which is just a generalized Lamb-shift 
Mi- In the closely related case of the interaction be¬ 
tween two driven Rydberg atoms, this phenomenon origi¬ 
nates the van der Waals blockade of the Rabi oscillations. 
This is the idea behind neutral atoms quantum gates 0, 
where the excitation of the target atom is blocked as the 
energy shift of its Rydberg level exceeds the value of the 
bare Rabi frequency. However, recent findings suggests 
that other dynamical effects might be relevant under cer¬ 
tain conditions-cf. Refs. jTTl. n~2l|. 

On the other hand, state of the art techniques on 
atomic interferometry have proved useful in measuring 
accurately the dynamical phase shifts accumulated by 
the wave function of coherently driven atoms |[l3]. Those 
phase shifts contain information about the interaction of 


the atoms with the environment. This is at the root of the 
proposals of Refs. PHI to measure the interaction of an 
alkali atom with a macroscopic surface at submicron dis¬ 
tances. At this distance, the dominant interactions are 
expected to be the CP interaction and the gravitational 
interaction. In the experimental setup of Ref. ^4] the 
atoms are trapped in a ID vertical lattice, which allows 
for an accurate control over the relative position of the 
atoms w.r.t. to the surface. The uniform Earth gravi¬ 
tational field creates a ladder of localized states referred 
to as Wannier-Stack states. At first glance, the net ef¬ 
fect of the Casimir interaction reduces to additive shifts 
on the atomic energy levels, both external (i.e., those in¬ 
volving the center of mass d.o.f.) and internal (i.e., those 
involving the electronic energy levels) [h|. As for the 
case of the van der Waals blockade, this approximation 
assumes implicitly that the Casimir energy can be inte¬ 
grated out a priori in the energy levels, before the atom is 
driven. However, dynamical effects may deviate from this 
assumption. Let us take as an example an atom driven 
through a lambda configuration [see Fig^a)] in the pres¬ 
ence of a material surface. On the one hand, the atom is 
pumped by two Raman lasers from two low-lying states, 

| g) and |e), to common excited intermediate states, |i). 
It is the combined action of both lasers that gives rise 
to the effective coherent oscillation of the atom between 
the states |g) and |e) fl7l-(l9i|. On the other hand, the CP 
interaction of the atom with the surface results in similar 
excitation/desexcitation transitions which are mediated 
by virtual photons rather than laser beams -see Fig|2K<i). 
Therefore, it is possible that both processes interfere with 
each other affecting the overall dynamics of the atom. It 
is our purpose to investigate the conditions under which 
the interplay between the Raman interaction and the CP 
interaction provides appreciable consequences on the co¬ 
herent dynamics of an atom. We will show that, in the 
perturbative regime, the net effect is an effective shift of 
the Rabi frequency. 

In this article we concentrate on scenarios similar to 
those of Refs.fl, [13, E3) where atoms are driven through 
a combination of Raman lasers and microwaves at the 
time they interact with a macroscopic surface. Our ap¬ 
proach is based on the time evolution of the atomic wave 
function. As long as the surface resonant frequency is far 
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from the atomic resonances, the classical treatment of 
the EM response of the surface is expected to be a good 
approximation. The case of the interaction between two 
driven Rydberg atoms is left for a separate publication. 
The article is organized as follows. In Sec. |TT] we review 
the essentials of both the Rabi model and the Casimir in¬ 
teraction, and motivate our work. In Sec. uni we address 
the renormalization of the energy levels and of the laser 
vertices, and show the role of the non-additive Casimir 
terms on the coherent dynamics of a driven atom. We 
give explicit expressions for the effective shift of the Rabi 
frequency. Sec. IIV Al contains an explicit calculation of 
this shift for a 87 Rb atom oscillating between two de¬ 
generate Zeeman sublevels close to a perfectly reflecting 
surface. In Sec. IIV Bl the calculation is done for a s 'Rb 
atom oscillating between two Rydberg states. In Sec. |V] 
we present our conclusions. 
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FIG. 1. (a) Schematic representation of a lambda system, 

where two Raman lasers of frequencies uj p and ui a couple the 
states \g) and |e) to a third intermediate state |i) with cou¬ 
pling strengths QF gi and £l s ei respectively. ( b ) Schematic rep¬ 
resentation of a ladder system. In this case the intermediate 
state |i) lays between | g) and |e). 


II. ESSENTIALS OF THE EM INTERACTIONS 


A. Laser fields interaction. Bare Q Rabi oscillations 

Let us consider first the interaction of the atom with 
the external electric fields of two monochromatic lasers. 
An atom in free space with atomic levels {|z)} is described 
by the free Hamiltonian given by 

H^=Y^mi a) 

i 

from which the unperturbed time-evolution operator 
reads, U at (t) = The Hamiltonian of in¬ 

teraction of the atom with the electric field of lasers of 
frequencies ui s and u) p is 

= hY^M p gi cos uj p t\i)(g\ +n s ei cos ui s t\i)(e\ 

i 

+ h n s gj cosU s t\j){g\ + h nrm cosu p t\m){g\ 

j^i,m m^i,j 

+ h.c. (2) 

where the strengths of the driven transitions are f2^? er = 
— (g, e|d|r) • Ep iS , with d the electric dipole moment op¬ 
erator and E PiS the amplitude of the electric field of the 
lasers p and s at the position of the atom. The interme¬ 
diate states labeled by i are commonly accessible from 
\g) and |e). On the contrary, states labeled by j are only 
accessible from | g) by the action of laser s and states 
labeled by m are only accessible from |e) by the action 
of laser p. Adjusting conveniently the detuning of the 
lasers w.r.t. the transtion frequencies to the common 
states, with = u>i g ,i e ~w PiS , as well as the strengths 

of the transitions, it is possible to make the atom os¬ 
cillate coherently between the states | g) and |e). That 
is for instance the case of an atom in either a lambda 
or a ladder system like those of Fig[2] Provided that 
|Agi s ei | |QgW| Vi, the population of the common in¬ 
termediate states can be eliminated adiabatically 0 and 
the effective dynamics of the atom reduces to that of 
a two-level system [3]. Straightforward application of 
time-dependent perturbation theory with the interaction 
Hamiltonian of Eq. J2| yields the effective Rabi Hamilto¬ 
nian Hr = H^ + Hr 1 , with 

Hr = h(io g + Suj g )\g)(g\ + h(uj e + 6w e ) |e)(e|, (3) 

H™ t (t)=Y^n i e i “ Lt \g)(e\ + h.c., (4) 

i 

where Su) g e are the frequency light-shifts, 

H = Ei(^i) 2 /4ASi + and 


We briefly review here the two EM interactions which 
govern the atomic dynamics. These are, the interaction 
with external electric fields and the interaction with the 
vacuum EM field. 


1 Throughout this article we refer as bare to all those observables 
which are computed in the absence of surface induced quantum 
fluctuations. 

2 The probability of excitation to the state |i) is proportional to 
(n p ./A p .) 2 sin 2 A p t + (n s ./A s .) 2 sin 2 A a . t. 

' gi< gx > gi v ex < ex' ex 
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. = £<TO74A-, + E„^(^ m ) 2 /4Af m re¬ 

spectively, flj is the effective bare Rabi frequency 
associated to the transition to the common state (*), 
^ = <>!;,<>',(a;;, + A£)/4A£AA, and WL = 
is the effective laser frequency. We note that for a 
2-photon transition in ladder-configuration we must 
take lo s < 0 in all the above equations. The diagrams 
contributing to the effective vertices of interaction and 
to the light-shifts are represented in Figsj2ja) and (b) 
respectively. In particular, the equation for the effective 
vertex represented by the lower diagram of FigJ2)(a) 
reads 

/ dt'e - *"*^ _t/ ) HI™*\ie(t')e~ iUet ' 
i 2ri J o 

X HIT+ [e o (5) 

where U(f) =T-exp/ 0 t dt’[H^ + H l e ^(t')\ and we have 
kept only the leading order terms in the second row. 
In Eq. (j4j), far off-resonant terms and terms of order 
A^W/wje.ig smaller have been discarded. 



FIG. 2. (a) Schematic representation of the effective vertices 
of the Rabi Hamiltonian, H^ t \ eg (t) and H l R t \ ge (t). (6) Dia- 
gramatic representation of the light-shifts, Soj g and Su} e . 

Under the action of the Rabi Hamiltonian, an atom ini¬ 
tially prepared in a linear combination of states | g) and 
|e) will undergo coherent Rabi oscillations. This problem 
is solved since long ago fTjj ]. for an extensive review see 
[20l ]. For the sake of comparison later-on we are inter¬ 
ested in two particular cases for which analytical solu¬ 
tions are well-known. These are, that of equal laser fre¬ 
quencies in A-condiguration, u>l = 0, and that for which 
the rotating-wave-approximation (RWA) is applicable @. 
The solution of the former is obtained by taking u>l = 0 
in the equations of the latter. 

Transforming Hr into the rotatory frame with rotation 
matrix R = —u)i\g)(g\ and discarding counter-rotating 


3 The RWA is valid as long as the effective laser frequency is close 
to the atomic transition, | ujR—uj eg | <y u) eg , and much larger than 
the inverse time of observation, ljr T —1 . 


terms, the eigenstates and eigenenergies of the trans¬ 
formed Hamiltonian, e lRt [HR — ?iR]e _iRt , read 


|+) = cos 0 c \g) + sin0 c |e), (6) 

|-) = sin6* c |g) - cos6» c |e), (7) 

E± = ( E e + Eg + Hujl)/2 ± — Hr, ( 8 ) 

A IHI 

with cos20 c = ——, sin20 c = ——, (9) 

Hr Ur 

A = OJL - Weg ~ 5ut e + Sujg, (10) 


n = h r = ^/a 2 + \h\ 2 . ( 11 ) 

i 

The wave function d/(T) for a time T > 0 of an atom 
driven by the above Hamiltonian and intially prepared 
at t = 0 in a linear superposition of the states | g) and 
|e), ’F(O) = a g (0)|5)+a e (0)|e), is given in the Schrodinger 
picture by 

vk(T) = U*(T)*(0) = [ 9 Uf (T)a s (0) + s Uf (T)a e (0)]| fl ) 

+ [ e Uf (T)a s (0) + e Uf(T)o e (0)]|e), 

where the components jU^(T), i,j = g, e, are given by 

g Xjf(T) = e-^~ A ^ T [co S (H r T/2) - t A sin (HrT/2)], 

Hr 

g Uf(T) = _j e -*K- A / 2 )rMsin(HflT/2), 

Hr 

e Uf(T) = e -*(“e+ A /2)Tjco S (HrT/2) + iA S m(n JJ T/2)], 

Hr 

eU f(T) = _ ie -<(<-.+ A / 2 )rM sin qi rT /2). (12) 

Ur 

B. Vacuum field interaction 

We evaluate now the effect of the vacuum fluctuations 
on the dynamics of a free atom for the case that the state 
of the atom is a coherent superposition of the states |<?) 
and \e) and the vacuum flucutations contain the inter¬ 
action of free photons with a closeby dielectric surface. 
In this respect, we use a semiclassical approach based on 
linear response theory. It consists of considering the pho¬ 
tonic states as dressed by the classical interaction of free 
photons with the dielectric surface. This implies that the 
linear response of the EM field i.e., its Green function, 
includes the scattering with the surfac^E The interac¬ 
tion of dressed photons with the atom is treated quan¬ 
tum mechanically. As we did above for the driven atom, 
we restrict the calculation of the time-evolution operator 
to the subspace {\g), |e}}. 


4 The resultant EM interaction is also referred to in the literature 
as ’body assisted’ interaction [jjj. 
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We constrain ourselves to the dipole approximation, 
W = — d-E(R) — m-B(R), where d and m are the atomic 
electric and magnetic dipole operators respectively and 
E(R) and B(R) are the electric and magnetic quantum 
field operators at the location of the atom, R. They can 
be decomposed in terms of oj modes as 

pOO 

E(R) = / dw[E(R;w) + E t (R;o;)], 

J 0 

poo 

B(R) = / dw[B(R;w) + B t (R;o;)], 

Jo 


induced by the transition dipole moments (g|d|e) and 
(< 7 |m|e) is negligible. We assume this condition in the 
following and apply time-dependent perturbation theory 
f2ll ] for the calculation of the time-evolution operator pro¬ 
jected on the subspace |<E>) = {|e), |g)} <g> | 0 ), 

U g(T) = [|6)(0| ® (\g)(g\ + |e)(e|)] U 0 (T) (15) 

x T-exp f dt Uj(t) W Uo(t) 

Jo 

x [(lfl)<ffl + |e>(e|)®| 0 )< 0 |]. 


where E(R; w)[B(R; w)] and Et(R;w)[B*(R;u;)] are the 
creation and annihilation electric (magnetic) field opera¬ 
tors of photons of energy huj at the position of the atom, 
R. The quadratic vacuum fluctuations of E(R;w) sat¬ 
isfy the fluctuation-dissipation theorem (FDT) relations 
at zero temperature ( 22 | . 

_2 

(0|E(R; lo ) ® Et(R / ; w)|0) =-=-9f[G(R, R'; w)l, 

7T6qC z 

(13) 

(0|B(R; w) ® B t (R / ; w)|0) = -^3[V fl AG(R,R';u) A 


1. CP interaction in the non-degerate case. Atomic level 
shifts 

In the non-degenerate case, with co eg T 1, at any 
order in T the diagrams which weight the most in the T- 
exponential of the diagonal components of and 

e U^, are those in which the atom transits through inter¬ 
mediate virtual states before getting back repeatedly to 
the original state, g and e respectively. Each transition 
through intermediate states is accompanied by the emis¬ 
sion and reabsorption of a single virtual photon which is 
reflected off the dielectric surface [see Fig[3]?j)]. 


where |0) is the EM vacuum state in the presence of the 
surface, 3 denotes the imaginary part and G(R, R';w) is 
the Green function of the Maxwell equation for the EM 
field, 




- V A (/i, 


-l 


Z, Z' > 0. 


V) A]G(R, R'; u) = d (3) (R, R')I, 

(14) 


In this equation e r and g r are the relative electric per¬ 
mittivity and magnetic permeability tensors respectively, 
and the atom’s position lies to the right of the surface, 
Z > 0 [FigEIKa)]. The Green function can be decomposed 
into a free-space component and a scattering component. 
The contribution of the free space term to the Casimir 
energy is the ordinary free-space Lamb-shift that we con¬ 
sider included in the bare values of the atomic transition 
frequencies. 

Next, we consider IF as a perturbation acting upon 
atomic and dressed photon states. The time-evolution 
operator of the latter is Uo(f) = [U at ® U 7 ](f), with 
U Ht) = Er=in“= 1 E 7 kj ^.e-^‘|7k„e,)(7 k 3 ,e i | the 

time-evolution operator of dressed multi-photon states of 
frequencies uij , effective momentum fi.k j and polarization 
Cj. In the following , we will refer to W as Casimir-Polder 
(CP) interaction. 

We will show that for the case that the doublet 
{| < 7 ), | e)} is non-degenerate in comparison to the obser¬ 
vation time, Lo eg T -1 , the net effect of the vacuum 
fluctuations is an atomic level shift . On the contrary, for 
uj eg T we will show that, in addition, the vacuum 
fluctuations may induce Rabi oscillations in the degener¬ 
ate doublet. In the latter case, the perturbative nature of 
the calculation is preserved as long as the CP interaction 




FIG. 3. (a) Atom (in red) placed at R in front of a plane 

surface perpendicular to the Z axis, with relative permit¬ 
tivity e r and permeability g r . (b) Feynman diagram of a 
multiple reflection process which contributes to 9 U^ in the 
non-degenerate case, u> eg T 1. The diagram equals three 
one-particle-irreducible additive phase shift factors, ( 5E gg ) 3 
-dissipative factors are omitted here for brevity, (c) Feyn¬ 
man diagram of a single reflection process which contributes 
to the non-additive shift factor, 5E ge . Left-handed (L) and 
right-handed (R) photons are created/annihilated at the time 
the atom transits from \g) to |i) and from |i) to \e) re¬ 
spectively. ( d ) Feynman diagram of a multiple reflection 
process which contributes to e U^ in the degenerate case, 
oj eg T ^ 1. The diagram equals three one-particle-irreducible 
non-additive phase shift factors, ( 5E gf .) 2 5Ef g -dissipative fac¬ 
tors are omitted here for brevity. 

On the contrary, the diagrams corresponding to the 
off-diagonal components, gU^ and e U^, weight much 
less and the T-exponential can be truncated at leading 
order. That involves a single diagram in which the virtual 
photons are created at the atom in one of the states and 
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annihilated at the atom in the other state [see Fig|3](c)]. 
Adding up the corresponding diagrams we obtain the 
components, 


g \jW (T) = e -iK+6-bE» s )T-r» s T/2 ] 

9 Uf (T) = —e~ iUgT [8E ge - ihT 9 ge /2]/E eg 
+ e~ iWeT 8E ge /E eg , 


e U Y{T) 
e U Y(T) 


= e -i(.u«+h- x 8Et e )T-Tl e T/2 _ 

= e-^ T [5El g - ihTl g /2]/E eg 
- e~^ gT 8E 9 g /E eg . 


From these expressions we deduce that the off-diagonal 
components of can be neglected up to terms of order 
0{8E 9 ’ gge /E eg ) for oj eg T 1. Thus, the effect of the 
CP interaction reduces here in good approximation to 
a renormalization of the energy levels of the atom. In 
particular, E 9>e -+ E g , e + 8E^ ee - ihT 9 g f e J2. 

In the above equations the additive energetic and dis¬ 
sipative terms have the usual expressions, 


SE S, 


T 9 

99 


- e • <16) 

*>7k,€ y 

% E ®(u g imAW\gM 2 S(u g i-u).(17) 


Analogous expressions hold for 8E £ e and r® e with the 
substitution e o g. 

The non-additive energetic terms are, 


= - E 

^>7k,€ 

(e, 0|W|i, 7)(z,7|W|<7,0) 

HU + flUJig 

SE % = - E 

{ 9,0 |W|i,7)(z,7|W|e,0) 

Hut + huJie 


( 18 ) 

(19) 


As for the non-additive dissipative terms, they are 

9-77- 

r s e = i 2 E 0 Ki)(ff 1 6|VF|7,i)(7,*|IF|e,O)(5(a;g J -ui), 

*,7k,e 


and an analogous expression holds for r® g with the sub¬ 
stitution e g. Single superscripts, g or e, in the ex¬ 
pressions for r and SE denote the reference frequency for 
the transitions within the sums. Double subscripts, gg , 
ee, eg or ge , denote the bra and ket states in the quan¬ 
tum amplitudes. In Appendix [A] we evaluate the electric 
and magnetic vacuum field fluctuations which enter all 
the above quantities via the FDT. 


2. CP interaction in the degenerate case. CP induced Rahi 
oscillations 

For uj eg T < 1 we review the calculation on Ref. 0 
adding up some details to it. In this case the diagrams 
which weight the most, both in the diagonal and in the 


off-diagonal components of are similar to those of 
FigEfb) i n the non-degenerate case, but for the fact that 
after each emission-reabsorpsion of a single virtual pho¬ 
ton the atom may arrive at either state, | g) or |e), with 
a similar probability. Fig[3][c) illustrates this process. 
Initially, a virtual photon of left-handed circular polar¬ 
ization, 7 L , is created at the time the atom is in state 
| g). Later, that photon is reflected off the surface turn¬ 
ing into righ-handed circularly polarized, r y R . Finally, the 
photon is annihilated at the time the atom gets to state 
|e). Diagrams involving n-photon intermediate states 
with n > 1 can be disregarded in good approximation 
since their contribution is of the order of (Tw) 1_n times 
smaller than the diagram with n single-photon interme¬ 
diate states, with u> the typical frequency of the vacuum 
photons. In the non-retarded regime, ui 1/T, so that 
multiphoton intermediate states are clearly negligible. In 
the retarded regime, u ~ c/R , so that their neglect is 
possible for T R/c, which is a realistic condition too. 
A typical diagram which contributes to the off-diagonal 
components of is depicted in FigEKd). Their sum¬ 
mation yields the following recurrent formulas, 


(T) = e~ iueT [l + E^IUf (T)], (20) 

n= 1 

e Uf (T) = e~ iui£T E (T), with 

n =1 

E) u rcn = ~ ,r ^ 1)! K^ - ihT ' g / 2)^ ( "E<(t) 

+ (SE £ ge -ihT £ ge /2)8^- 1 l\]^(T)], 

§ (n) v w (T) = ~ lT ^~ 1)! [(^f e - ihTt/2)5^- l l\]^ (T) 

+ (^fg-^rfg/2)^”- 1 )uf(T)], 

S {1 IK(T) = - iTh ~ X (* E 9 9 - ihT£ 9 9 l 2 )> 

6 {1 IK ( T ) = -iTh~\5E £ eg - ihT%/2\ ( 21 ) 


where we have used £ = ( E g + E e )/ 2, uj£ = £/h. Analo¬ 
gous expressions hold for eU^ (T) and (T) exchang¬ 
ing the subscripts ef>jin the above equations. 

From the diagram of Fig(3]7i) we observe that each fac¬ 
tor {8Ef g — ihTf g /2 ) flips the state of the atom from |e) 
to | g), while each transposed factor produces an opposite 
flip. This is analogous to the action of the factors Ml/2 
and Ml* /2 of the Rabi Hamiltonian respectively, except 
for the fact that here the damping terms break the time 
reversal symmetry. As a result, possesses the same 
functional form as in Eq. m, with the following sub¬ 
stitution of the bare parameters in Eq. m with the tilded 
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ones defined below @, 

^9 = + SEgg/h — iTgg/2 , 

We=W e + 6E% e /h- *rf e /2, A EWe-W 9 , 

« = 2 $Eg e /H - *rf e , n* = 2 <S£*ft - *rf g , 

|q| 2 eeqq*, = y^Qp + A 2 . 

This means that the Casimir-Polder interaction may 
indeed induce Rabi oscillations between two quasi¬ 
degenerate states Ea¬ 
rn. VACUUM-INDUCED SHIFT ON THE RABI 
FREQUENCY 

In this section we derive an expression for the total 
shift induced on the Rabi frequency of a driven atom 
by vacuum fluctuations, SHr = 6 iHr + S 2 HR + S 3 QR. 
In this equation we distinguish three different kinds of 
shifts, namely, that induced by the additive CP terms, 
SiHr-, the one induced by the non-additive CP terms, 
&2 Hr; and the one induced by the renormalization of the 
vertices of interaction in 63 Hr. 


A. Shift by additive CP terms 


It is obvious that the additive terms of the Casimir in¬ 
teraction cause a shift on the atomic levels which induces 
a variation on the Rabi frequency [24j]. First, the effec¬ 
tive couplings Hi appearing in Eq. Q experience a varia¬ 
tion S a Hi as a consequence of the shifts of the detunings. 
Second, the global detuning of Eci. flTUl) . A, changes by 
an amount 5 A due to the shifts of the levels E e and E g 
as well as to the variations induced on the light-shifts, 
Soj e and Su g [see Fig|3Ja)]. As a result, as a function of 
the bare quantities f2, Hr and A, and of the additive CP 
terms, we have at leading order 


SiHr ~ —— y 5 a tti -I- — 5A 

Hr Hr 


-n 


\W P giKi , SE u ~ SE$g 8EI-8E* 


4hH R A % 


A s2 


HPgl(SEi - 8E° gg ) {6EI - 6E&). 


+ 


4MIr 

A 


A p2 


A s2 


ftQ 


R L 


H s2 

go 


SE’ g -SE’„- Y -^(SE^-SEj,) 

j^i,m 


93 


C)p2 

E T 


m^i,j 


4A p2 


( 22 ) 


5 Note here that Cl* so defined is not the complex conjugate of £7 
because of the damping factors. 


B. Shift by non-additive CP terms 


We have found in Sec. Ill B 21 that for uo eg T < 1 an in¬ 
duced Rabi frequency may be provided by vacuum fluctu¬ 
ations. In the case of a driven atom, |Q| ~ 1/T, so that 
the degenerate condition is equivalent to the so-called 
deep strong coupling (DSC) regime, |f2| u} eg [H]. The 
corresponding shift on H is given by diagrams in which 
pairs of consecutive Raman vertices and pairs of consec¬ 
utive CP vertices alternate -eg. diagram of FigQJfr). As 
an example, we compute the leading order terms which 
contribute to the variation of S g U 9 . To this aim, 

we treat W + Hr* as a perturbative interaction upon 
Hq*. Leading order contributions are of cubic order and 
correspond to the terms H’^W 2 and W 2 Hr*. At this 
order, we find 


W ~ — {iT/2E eg ) 


—iu) e T o* 


H*(SE e qe - 


° T H{5Etg 

iW e ge /2) 


ihTlg/2) (23) 
for uj eg » |Q|, 


^ ■ ihr9 ^ /2) (24) 
+ fi* (SE ge - *ftr ge /2)] for Weg < |fi|, 


where, in order to simplify matters, we consider zero 
global detuning, A ~ 0, and we assume that all the tran¬ 
sition frequencies to intermediate states are much larger 
than |f2|. Straightforward comparison with the equation 
for g U^ reveals that Ea. (l24l) is its term of 0(T 2 ) for 

n->n + 2(«5f e /ft-*r® e /2), 

H* ^H* + 2{5El g /h-iVl g /2). (25) 


A similar calculation at order four in W+Hr* contain¬ 
ing the terms {H^fW 2 , W 2 ^*) 2 and H^W 2 , 
yields the following result, 

W = e-^ Tl -^fl(8E° g - 3hTlJ2) 

+ {8E e ee - ihT e ee /2)/2], (26) 


which holds for |f2| either much larger or much smaller 
than w eg . Again, comparison with the equation for g U^ 
reveals that Ea. tE^l) is its term of 0(T 3 ) for 


A -a A + [{5El g - 5E e ee )/h - i(r® fl - hT% e )/ 2 }. (27) 

As expected, additive terms appear always as energy 
shifts irrespective of the magnitude of the ratio \H\/ui eg . 

Disregarding dissipative terms, the shift of A was al¬ 
ready accounted for in Eci. (1221) . The additional CP in¬ 
duced shift on the Rabi frequency is therefore 

2Q 

^Hr ~ A— 9?{AE eg },for u) eg < |Q|. (28) 

“A 

Alternatively, the shifts on the Rabi frequency induced 
by the additive and non-additive CP terms SE gg , SE^ e , 
SE ge , SE^g for uj eg <C |Q| can be computed out of the 
renormalization of the eigenstates |+), |—) and their cor¬ 
responding eigenenergies -see AppendijBl 
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FIG. 4. (a) Diagrammatic representation of a multiple reflec¬ 
tion process in which additive CP terms induce a shift on the 
Rabi frequency of the kind of Ea. (l22l) . One-loop-shifts 

in the levels E g , E e and Ei are included. ( b ) Diagrammatic 
representation of a multiple reflection process in which non¬ 
additive CP terms induce a shift on the Rabi frequency of the 
kind 52^r of Eg. (1281) for uj eg |fi|. (c) Diagrammatic rep¬ 
resentation of the one-loop vertex renormalization processes 
which generates SHf g and 5Hf g respectively. They give rise 
to a shift on the Rabi frequency of the kind £ 3 Qr of Eg. (1331) . 
The transient states i, g and e belong to the same energy 
levels as the states i, g and e respectively. 


C. Shift by vertex renormalization 


with 


= E 


-SI?,, 


91 


2e 0 c 2 l 


‘■'isTrjdig • G (u>i g ) ■ d ?; J 


tgJ u igi 


+ - 


d u 1 r 


n J n u 2 + uSf, 




■Tr{dig • G (iu) • d 5 } 

TU, is 

j f°° du U 2 (U 2 - UJjgUJei) 

ne 0 c 2 J 0 (u 2 + u 2 g ) (u 2 + w 2 0 ) 


x Tr{djg • G (iu) ■ d-J, 

where d ah denotes (a|d|6). An analogous equation for 
SH™*] ie (t) yields 

SH™* lie (t) — COS UJpt + SH- e cos u s t, (30) 


with 


S H U = E 2^ [“eiTr{d iS ■ G*(w e j) • d-J 


+ 


d u u 2 
u 2 + w 2 


= E 


-sir. 


CqC L U>i g - W e 


-Tr{d,e • G(m) • d-J 


-Tr{d,g • G*(w is ) • d-J 


w; 


+ 7T 


'L 


-Trjdjg • G (uei) ■ d-J 

du U 2 {U 2 - UJigUJei) 

0 w+l 


x Tr{dig • G (iu) ■ d-J 


Lastly, we are left with the effect of vacuum fluctua¬ 
tions on the Raman vertices appearing in Eq.Q. This 
corresponds to diagrams in which single Raman vertices 
and single CP vertices alternate. At leading order of 
time-dependent perturbation theory, the variations cor¬ 
respond to the diagrams of Fig(4)(c) . They read 

6Hlf\ ig (t) = (t,6|U“ t (i)U w +(i)flS t (*)U w (t)U ott (t)|fl,6) 

-w™m) 

-n- 2 V f dfV( T -*'X“ +w s)e it,w J*,6|lF|s,7) 

~.J 0 

7tj J 

Jo 

+ /r 2 V [ dt , e i ^ T ~ t '^ u+Ue ^e it uli (i, 0 |W|e, 7 ) 
~.Jo 

xe-^-Je, 7 |fl£*(i)|i, 7 ) 

rt 

x / (1, 7 |W|< 7 ,0) 

Jo 

~ cos w J + (5iL*g cos u a t, (29) 


In the above equations the tilded states |z), |g) and |e) 
belong to the same energy levels as the states i. g and e 
respectively. For simplicity, we have assumed equal en¬ 
ergies for all the intermediate states, w- = Wj Vz. Far 
off-resonant terms w.r.t. to the transition \g) o \e) and 
rapidly evanescent terms have been discarded. The com¬ 
plete expressions for 8H™ t \ ig (t) and Affgjjjji) can be 
found in Appendix [C] We note that all the terms above 
are of the order of ^ g s eg 8E^’ g e ee / 0Ji g ,ei- Therefore, their 
contribution is generally relevant for very strong CP in¬ 
teraction w.r.t. the transition frequencies to intermedi¬ 
ate states. This might be for instance the case of an 
atom which is made oscillate between two close Rydberg 
states in ladder-configuration, near a metallic surface. 
Thus, ui e i > 0 has been assumed in the above equations. 
It is now straightforward to calculate the effect of ver¬ 
tex renormalization (n) on the coherent evolution of the 
atomic wave function. By comparing Eas. (129l30D with 
Eq.([2]) we read that it consists of a shift on the effective 
bare Rabi frequencies. Generally we find, 

SvSU = (6H p gi /hCl p gi + 8H s ei /m s ei )Q,i, (31) 
while for the special case |2w PjS — wj -C |A| we have, 
SvOi = m p gi + SH s gi )/hn p gi + ( 6H° ei + 

(32) 


































It is now straightforward to compute the total variation 

of the Rabi frequency due to vertex renormalization, 4 2 d 3 / 2 




R — 


n 

D R 


y s v fi z . 


(33) 


IV. CP INDUCED SHIFTS ON THE RABI 
FREQUENCY OF A RUBIDIUM ATOM CLOSE 
TO A REFLECTING SURFACE 


A. Oscillations between degenerate Zeeman 
sublevels in A-conflguration 


In Ref. Em it was found that a CP induced Rabi fre¬ 
quency as large as 27 tHz could be obtained between the 
Zeeman sublevels | g) = \h 2 S\/ 2 ,F = 1,to_f = —1) and 
|e) = \5 2 Si/2,F = l,m_p = +1) of a s 'Rb atom at zero 
temperature when it is placed in the vicinity of a perfectly 
reflecting surface parallel to the quantization axis, and 
as long as | g) and |e) remain degenerate, u eg -C 27tHz. 
Symmetry considerations imply that the CP interaction 
respects this condition since the additive energy shifts 
of both states are equivalent. However, the presence 
of stray magnetic fields may induce a Zeeman splitting, 
Huj eg = I^-bBo, with B 0 the effective strength of the ma- 
gentic fields along the quantization axis, say X, which 
could break the quasi-degeneracy of the states and vio¬ 
late the condition ou eg <C Qr = 2SE^ g /h. In order to 
avoid this potential problem we propose an alternative 
setup in which we aim to measure the CP induced shift 
on the bare Rabi frequency of a driven atom in the DSC 
regime. The computation is in all points equivalent to the 
one of Ref. & The advantage of working with a driven 
atom is that, for a sufficiently high value of the bare Rabi 
frequency, |fi|, the unknown quantity ui eg contributes to 
the actual Rabi frequency with a shift ~ uj 2 g /2\VL\ 1 which 
can be made negligible for sufficiently large |Q| w.r.t. the 
shift provided by the non-additive CP terms, 2 SE^ g /h. 

The setup is sketched in Fig|5l where the parameters 
have been chosen so that the value of the effective bare 
Rabi frequency is Qr = |Q| = 27 t 20 Hz. Two Raman 
lasers of equal frequency and opposite circular polariza¬ 
tion drive the transitions from |g) and |e) respectively to 
a virtual state close to |ii) = |5 2 P 1 /2,P’ = l,mp = 0) 
and |i 2 ) = |5 2 Pi/ 2,P = 2,m_F = 0). For convenience 
we take A s gii = A p eii = A 0 = 10GHz, and the same 
intensity for both lasers, |Qo| = 27 t 8.0 MHz, such that 
|Oo| <C A 0 <C i^Is5P 1/2 - Denoting their electric field 
strength by | E 0 |, we define 

Qo = fi- 1 (5 2 5 1/2 ||d||5 2 P 1/2 )|Eo|, (34) 

and the effective bare Rabi frequency can be computed 


5 2 P3/2 
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FIG. 5. Schematic representation of the action of two Raman 
lasers which drive a 8 ' Rb atom between two Zeeman sub- 
levels of the ground state, \g) = |5 2 Si/2,F = 1, wf = — 1) 
and |e) = \5 2 Si/2,F = l,mF = +1). The quantization 
axis is taken along the X direction. The lasers have oppo¬ 
site polarization, a + and er _ , equal frequency, A ~ 795nm 
such that Ao = 27rl0.0GHz and equal intensity such that 
Do = 27 t8.0MHz. They couple both states to a pair of 
common intermediate states of the level |5 2 Pi/ 2), |ii) = 
|5 2 Pi/2,P = 1 ,m F = 0) and |i 2 ) = |5 2 Pi/ 2 ,P = 2, m F = 0), 
with transition frequencies u>l S5Pl/2 and to^g 5Pl/2 + 5u> 1 re¬ 
spectively, where u>l S5Pl/2 = 27r377.107THz and the hyper- 
fine interval is 6u>i = 27r0.817GHz. Grey lines represent 
secondary couplings of the three Zeeman sublelevels of the 
ground state, |5 2 Si/ 2 ,P = 1), to other states in |5 2 Pi/ 2 ). 
They are relevant for the computation of light-shifts. In 
addition to the D1 transition line, relevant to the calcula¬ 
tion of S2^1r are the D2 line transitions from | g) and |e) 
to the Zeeman sublevels of 5 P 3/2 with m P = 0. For those 
transtions, uj\s 5 p 3/2 = 27r384.230THz, Si 02 a = —27r0.072GHz, 
5id2b = 27r0.157GHz. Also, relevant to Si ftp are the transi¬ 
tions 5P!/2 65*3^2 and 5 P 3/2 4 D 3 / 2 . 


as a function of Qqi Aq and 5lo\ as, 


n = \E 0 \ 2 [(g\d + \i 1 }(e\d-\i 1 )*/2A 0 
+ (g\d + \i 2 ) (e\d-\i 2 )*/2,{A 0 + <5wi) 
|D 0 | 2 


24Ag 


Suj 1 , 


(35) 


where the hyperfine interval satisfies Suii <C Aq and the 
dipole moment operator is expressed in the spherical ba- 
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sis. The expectation values and transition frequencies 
have been taken from Ref. [25|. 

Beside the couplings between the states | g) and |e) and 
the common intermediate states |*i) and |i 2 ), the Raman 
lasers also couple the Zeeman sublevels of the multiplet 
|5 2 5 'i/2> F = 1) to all the hyperfine Zeeman sublevels 
within |5 2 Pi/ 2 ). The pairs of couplings are depicted in 
Figj5] with gray lines. The result is that while the dif¬ 
ferential light-shift between |e) and | g) is exactly zero, 
there exists a differential light-shift between the state 
| 5 2 Si/ 2 ,F = 1 ,771 p = 0) and the other two states. This 
light-shift breaks the degeneracy between the three Zee- 
man subleves and rises |5 2 <S' 1 / 2 ^ P = 1, tuf = 0) above | g) 
and |e) by an amount = 10Hz27rP In turn, this 

is an advantage, since the states | g) and |e) become the 
lowest energy states and no dissipative decay terms enter 
in the calculation. Only at very short distances the pos¬ 
itive magnetic energy shifts, ^/192neoc 2 Z 3 , may take 
the states \g) and |e) over |5 2 Si/ 2 ,P = 1,to_f = 0) a 

We proceed to compute <5 VLr. In the first place, 
given that the bare global detuning is approximately 
zero, the Rabi shift induced by additive CP terms is 
caused by the variation of Ao in the equation for Q 
[Eg. (1551) ]. The variation of A 0 is itself due to the dif¬ 
ferential level shift between the states 5 Si/ 2 and 5 Pl/ 2 . 
In the near field, Z < lOOnm, the main contribution to 
it comes from the difference between the reduced dipole 
matrix elements (6 2 5 , 1/2 ||d||5 2 P 1 /2), (4 2 P 3 /2||d||5 2 P 1 / 2 ) 
and (5 2 S'i/2||d||5 2 Pi/2 i 3/2) [H|. As a result, we have 
~ — SE 9 g ), where i stands for any of 

the intermeciiate states, i\ or i 2 . Using the equations in 
Appendix 0 for the additive CP terms as functions of 
the Green’s tensor which is, for a perfectly conducting 
reflector, 


'2ikZ 

G xx (Z\u>) = G yy (Z\u) = 39 ^ 2^3 (~1 + 2 ikZ + 4 k~Z 2 ), 

2ikZ 

G zz {Z ; w) = 167r £, 2 £3 ( _1 + 2ikZ), (36) 

the differential level shift in the near field can be written 
as 


5 El - 5 EL ~ 


-1 


l(6 2 *S' 1/2 ||d||5 2 P 1/ 2>| 2 


99 487re 0 Z 3 L 1 

+ E I ( 77 - 2 -P ) 3/21 |d| |5 2 P\/ 2 ) | 2 


(37) 


n—4 
-2 


- |(5 2 S 1/2 ||d||5 2 P 3/2 )| 2 - 2|(5 2 5 1/2 ||d||5 2 P 1/2 )| 


As for the Rabi shift induced by non-additive CP 
terms, the calculation was already carried out in Ref. 0- 
Here we just give the final formula, 


— 


1 


16Sc7r 2 e 0 -Z' 2 


dw- 


f(u) 


3/2J+3/2 

E 

j—l/2,F—\j—3/2\ 


F 

W 5S5P,- 


x (e\d- |5 2 Pj, F,rriF = 0)(g|c?+|5 2 P ) -,P,tof = 0), 



FIG. 6. Graphical representation of the Rabi frequency shifts 
as a function of the distance to the reflecting surface. 


where f(u) = e^ 2 “(l + 2u — 4u 2 ) and Kjp = ^issPjZ /c, 
with u}^ S5P . being the transition frequency from the hy¬ 
perfine level P of 5 2 Pj to the states | g), \e). 

In Fig [6] we represent the values of the Rabi frequency 
shifts as a function of the distance to the surface in the 
non-retarded regime, Z < lOOnm. It is clear that <5 2 Uf 
dominates over in our setup. 


B. Oscillations between Rydberg states in 
ladder-configuration 

Let us consider now the Rb atom close to the reflect¬ 
ing surface oscillating between two Rydberg states, | g) = 
|48 2 S’ 1/2 ,P = 1, rrij = -1/2) and \e) = |49 2 S’ 1/2 ,P = 
1, rrij = —1/2), under the action of a 7r-polarized mi¬ 
crowave source which is resonant with the two-photon 
transition, wo = w eg /2 = 2-7T 35.2387GHz [FigJT]. This 
setup is meant to mimic the situation of some hybrid 
quantum systems in which the coherent manipulation of 
Rydberg atoms is used to probe the EM fluctuations near 
a surface -eg. Ref.@|. The microwave source couples 
the states | g) and |e) to the intermediate states |H) = 
|48 2 P 1/2 ,m J = -1/2) and \i 2 ) = |48 2 P 3/2 ,rrij = -1/2), 
with equivalent detunings Ai = A ilSje = — 2 f 1.093GHz, 
A 2 = A i 2g}e = —27r 0.1611GHz, and with nearly equiv¬ 
alent strengths —^ il 2e . We set O 0 = 
2 27r 53.884MHz so that Qr = 27r 10.0MHz. We note that, 
despite the two-photon resonance condition, there ex¬ 
ists a bare global detuning due to the light-shift, A = 
Sojg — Suie = —2n 1.719MHz. 

As already done in the previous section, the com- 
putaion of S involves the evaluation of differential 
light-shifts. In the non-retarded regime this implies the 
addition and substraction of the squares of the reduced 
dipole matrix elements. While the nS states couple only 
to nP and (n — 1)P states in good approximation, nP 
states couple also to nD, (n — 1)D and (n — 2 )D states. 
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FIG. 7. Schematic representation of the action of a 7r- 
polarized microwave source which drives the oscillations of a 
87 Rb atom between two Rydberg levels, | g) — |48 2 5'i/ 2, w.y = 

— 1/2) and |e) = |49 2 Si / 2 ,Trij = — 1 /2). The quantization 
axis is taken along the X direction. The intermediate states 
are |ii) = |48 2 Pi/ 2 ,mj = -1/2) and |i 2 ) = |48 2 P 3 / 2 ,m ;/ = 

— 1/2). Exact a priori resonant condition is assumed, u>o = 
w eg /2 = 27t 35.2387GHz. Hence, A 9 ; 12 = A e q 2 . 


processes mediated either by two 7 r virtual transtions or 
by two consecutive er + and er_ transitions yield nonvan¬ 
ishing contribution. In Figs|8]A) and (c) we depict two 
of these processes. It turns out that the contribution of 
those processes involving 7r virtual transitions vanishes 
for Ei r ~ Ei 2 . Therefore, we are left only with processes 
of the kind of Fig HD c), 

_on 

(4 °> 

x [(48 Pi/2i +l/2|d+|g) + (48P 3 /2, +l/2|d+| g)\. 




This implies that the order of magnitude of the differen¬ 
tial shifts in Eg. (1771) is ~ AE|fp 48P . Making use of the 
data in Refs. 1127, |28| we obtain, 


£>i Qr — 


A 8.147 1 


r n 18.148 


Q.rZ 3 (pl m) fl R Z 3 {fj,m). 


2tt MHz 


18.282 
Z 3 (fi m) 


2tt MHz. 


(38) 


As for the shift induced by the renormalization of 
the laser vertices, S 3 CIR, we use Eg. (1771) combined with 
Eg. (1321) since there exists only one microwave source. 
Nonetheless, it can be verified that, in the non-retarded 
regime, SH^ 2g and 2e are much smaller than SH^ 2g 
and 6 H? , and hence negligible. This means that pro¬ 

cesses like that depicted in FigHDa) in which two vir¬ 
tual transitions take place between | g) and |e) can be 
neglected, and we end up with S v Cli 12 ~ — [SHf i + 
5H ? ]f2,/f2o- Discarding next the resonant components 

of SH ? 2g and SH^ 2e , which are much smaller than the 
non-resonant components in the non-retarded regime, 
and using the Green’s function of Eg. (1771) in the near 
field, we obtain 


SH ’“ " ElSWflfi-W 

- (3/2)(i\d-\g)(i\d+\g)\. (39) 

An analogous expression holds for SH^ e . In Eo. dTTl) con¬ 
servation of total angular momentum implies that only 


FIG. 8. Schematic representation of three processes which 
contribute to 8 H gi2 , (a); 8 H° i2 , (b); and SH p gii , (c), respec¬ 
tively. Wavy lines depict virtual transitions driven by vacuum 
photons. Dashed lines depict actual transitions driven by the 
microwave source. I 11 the non-retarded regime (a) is negligible 
and, after considering all kind of transitions, only processes 
of the kind of (c) survive. 


Finally, adding the contribution of SH? e and making 
use of the fact that the reduced dipole matrix elements 
between Si / 2 and Pj states hardly depend on j [23], we 
can write in closed form, 


83 ^ R — 
+ 


—17 217^ T 172 2 

f7p wo5127reo/LZ' 3 

|(48P||d||49S')| 2 ] : 


[|(48P||d||48S')| 2 
0.025 


Z 3 (^n 1 ) 


2t r MHz, 


(41) 


which is almost three orders of magnitude smaller than 
Jiflfl. A graphical representation is given in Fig[9] The 
ratio Sifla/Ssflii can be worked out from their expres¬ 
sions in Section eh We have, <5if7p/(i 3 f7p ~ A 2 /Auj 0 , 
with A a numerical prefactor of order unity. We conclude 
that, altough <5 3 f7p is generally a few orders of magnitude 
smaller than the ordinary <5if2p, it may have an effect on 
high precision measurements. 

We finalize this section with a comment on the setup of 
Ref. [l3l] . where a ®'Rb atom is made oscillate in ladder- 
configuration between two hyperfine-structure states, 

I g) = |5 2 Si/ 2 , F = 1 ,m F = 0) and |e) = \5 2 S 1/2 ,F = 
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FIG. 9. Graphical representation of the Rabi frequency shifts 
as a function of the distance to the reflecting surface. 


2, mF = 0), with io eg ~ 7GHz. The two states are con¬ 
nected by an Ml transition. An analogous calculation to 
the one performed above yields ~ U/ 100 Z 3 (nm), 

which is negligible for operational distances larger than 
lOOnm. 


V. CONCLUSIONS 

We have analysed all one-loop radiative corrections 
which contribute to the shift on the Rabi frequency of a 
driven atom close to a material surface. In addition to the 
shift induced by the ordinary additive variations on the 
atomic levels, SiEIr, two novel contributions have been 
reported. A shift induced by the non-additive Casimir- 
Polder terms, is found to dominte when the atom 

is made oscillate between two degenerate Zeeman sub- 
levels in lambda-configuration. A shift induced by the 
renormalization of the laser vertices of interaction, S^EIr, 
contributes at higher order than 5±EIr for an atom which 
is made oscillate between two Rydberg states in ladder- 
configuration. 
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Appendix A: Additive and non-additive CP terms 

In this Appendix we compile the expressions for the 
energy shift and dissipative CP terms, SE and T respec¬ 
tively, both additive and non-additive. As in Sec. IIIB1 
the single superscripts in the expressions T kl and SE 3 kl 


denotes the reference frequency for the transitions, j = 
g,e,E; while the double subscript denotes the bra and 
ket states in the quantum amplitudes, kl = gg , ee, ge, eg. 
In addition, we use the notation d a b = (a|d|6), m a b = 
(a|m| 6 ). We apply the FDT outlined in Section fllBI for 
the evaluation of vacuum field fluctuations at zero tem¬ 
perature. We separate electric and magnetic field contri¬ 
butions and, for the sake of simplicity, we assume that 
the surface possesses no chiral response, 


SEfa = 


53 [1 Ssj(Sie + <5jg)] 

i,Jk,e 


r(fc,6|Wej|i,7)(i,7|Wei|i,6) 


hui + hujij 


| (fc, 0|IF m |^,7)(^,7|I ; F m |Z,0) ^ 

I^u; + Hujij 

1 r°° 

= -2 F / dw w 2 531 1 - Ssj{5ie + Sig )] 

7T£0 C 2 Jo ^ 

x T ^ dfci • S[G(R, R; cc)] • dg 

(xJ + LUij 

-2 P / dw53[ 1 - 6ej{6ie + 6ig)] 

7re 0 C Jo . 

x Ti m fci ■S[V fl AC(R,R;ij)AVji]- m,j 

U) + Cdij 

27T 

^kl = ^2 55 + Sig)]Q(uJji ) 

*, 7k,6 


(Al) 


X 

+ 


(fc, 0|Wei|7, i)( 7 , i\W e i\l, 6 )6(uji - w) 
(k,0\W m \'Y,i){'y,i\W m \l,0)5(u j i - u 

2 53[1 ~ $Sj( 1 ~ — ^is)]0( w ii) tt ’jiTr{dfcj 

i 

3?[G(R, R;wji)] • d i; } 

^ c2 53[1 - Ssj(Si e + S ig )]Q(u}ji)TT{m. k i 
i 

3f[V A G(R, R; ojji) A V] • m^}. (A2) 


The factor [1 — Ssj(Si e +Si g )] in these expressions accounts 
for the removal of the states g and e from the sums when 
the reference energy level is £ = (E g +E e )/2 in the quasi¬ 
degenerate case. This ensures the perturbative nature of 
the calculation. 

In general, in the energy shift terms we can distinguish 
resonant (r) and off-resonant (or) components j 2 , 0 ] ■ 

The resonant components account for the single poles of 
the integrand in Eq. (ED, 


^kl\r ~ 2 531 1 fisjO- 

l 

x Tlr{d w -R[G(R,R;a; ii )]-d i ,} 

-^ 531 1 “ hj{&ie + <5i g )]0(wji)Tr{m fci 

eoc 

i 

■ 3ff[V A G(R, R; ujji) A V] • m^}. 


(A3) 
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For the off-resonant components, making use of the 
properties of the Green’s functions, G(R, R; — w) = 
G*(R,R;w), w 2 G(R,R; w) —> 0 as |w| —» oo, and em¬ 
ploying standard integration techniques in the complex 
plane [2| we find, 


5EI \°r = / du u 2 Y ,[1 - hi (*« + 8 ig )] 

Jo i 

d fci • G(R, R; iu) • d it 

x WjjTr- 5 —- 5 - 

u 2 + W 2 - 

d u Y,[l-h j (S i e + 8 ig )] (A4) 

m fci • SfV/j A G(R, R; iu) A Vr] • mu 
x w;,Tr-—=-. 


Finally we note that these expressions can be formally 
rewritten as functions of the atomic polarizabilities using 
the appropriate definitions Q. 


Appendix B: Renormalization of eigenenergies and 
eigenstates for oj eg -C |fl|. 


We consider the eigenstates of the Hamiltonian Hr 
of Eos. (13141) . |+), |—), as stationary states upon which 
W acts as a stationary perturbation. This is a good 
approximation for small effective laser frequency, wrT <C 
1, and for the case that the virtual transition between \g) 
and |e) are irrelevant in the CP interaction. Application 
of time-independent perturbation theory at order W 2 , up 
to 0(5EA 2 /\Q\ 2 ), yields the energy shifts 


SE = V 1 ( A 7 |^ I +)| 2 V l (- 7 | ty |+)| 2 

hck + Ei — E + ' hck — H\Q\ 

i^g,e,7k, e 7k,€ 

- - 2 (SE £ gg + SE £ ee ) + nSE%} 

+ m {SE ’’- 5EL) -m nsE -’ ) ' (B1) 

SE = sr Kb 7 l ^|-)| 2 ^ |(+, 7 n -)| 2 

, hck + Ei — E- hck + h\Q\ 

7k, e 7k,e 

~ l -{5E £ gg +5E £ ee )-K{5E £ eg } 

~ 2]f T\( 5E ° 9 ~ SE ^ + ( B2 ) 

As already explained in Section IIIBI the double sub¬ 
scripts in the quantities SE, gg , ee, eg or ge , denote the 
bra and ket states in the quantum amplitudes, while the 
superscript £ denotes the common reference frequency, 
cog = (ojg + w e )/2, for the intermediate atomic transi¬ 
tions involved in their calculations. We have assumed 
<C cui e ,LUi g ,LOi£ V* relevant in the sums over inter¬ 
mediate states, i ^ g,e, so that we can approximate 
SE £ gggee ~ SE £ eg }^)H. As anticipated in Sec. IIIIB1 


the net result of these energy shifts is a renormalization 
of the bare parameters which enter U^, 


OJg ^ UJg — LUg + SE gg /H, LO e ^ U) e — 

w e + 5E £ e /h, 

A-> A = A - (SE s ee - SE gg )/h, 

(B3) 

|nH|n| = |n| + 2R{5£f fl }, 

(B4) 

Gr -!t kin = + ^^{SE^} 

(B5) 

+ ^E £ gg -SE £ ee )-^{8E £ eg }. 



For the sake of completeness we compute the variations 
on due to the interaction W, <5U fl . To this aim we 
calculate the wave function at time T > 0 for the initial 
condition 'F(O) = | g), 

| $(T)) = [V R + SU RW ]{T)\g) (B6) 

= e -^ T (E++SE + ) j cos 2 ^ + gin 2(j c \e)/2] 

+ e~ lT( ' E ~ +SE ~' > [sin 2 9 c \g) — sin20 c |e)/2], 

where the renormalized (tilded) trigonometric functions 
are given by the expressions in Eq.© but for the re¬ 
placement of the bare parameters with the renormalized 
ones of Eas. (IB3HB51) . In turn, [U^ + 5U fl ] (T) presents the 
same functional form as the operator HJ fl (T) in Eq. m 
with the replacement of the bare parameters by the 
renormalized ones. For the sake of simplicity we choose 
A = —w eg , so that all exponential prefactors in front 
of the components of in Eq. m become e lul£T 

with uj£ = (ojg+LOe)/ 2. We obtain, at leading order in the 
energy shifts up to terms of the order of 0(ASE/h\£l\ 2 ), 

5 9 Uf(T)^e-^ T {cos(|f2|r/2) (B7) 

5E £ 

- i—f^[ cos (\n\T/2)T/2 + l^l” 1 sin (|n|T/2)] 

- ^{6E £ ge }sin(\Q\T/2) 

- ^[cos (\fl\T/2)T/2 - Ifip 1 sin (|n|T/2)]}, 

5 e \J R (T) ~ -e~ iU£T [i sin (|fi|T/2) (B8) 

+ (6Ei+6E £ gg )sm(\n\T/2)T/2h 

+ ^{5E £ g }T cos (\Q\T/2)}. 


Appendix C: Vertex renormalization 

We give the complete expressions for the one-loop 
vertex shifts, 5H^\i g {t) an d 8H l ^\i e {t). We restrict 
ourselves to the electric dipole approximation and as¬ 
sume that the atom is driven in ladder-configuration, 
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^iai ei ^ 0- 


SH^Ut) = E 4^ [wfeTr{d ie- • G*(w ie ) • djje^* 
i,e 

+ w ie Tr{d ie - • G(w ie ) • 


-fi p 

= EttE ^Tr{d l§ . G*( Wls ) • dije^* 


4e 0 c 2 L 

+w is Tr{d i§ • G(w ifl ) • d io }e < <"’- 2 “ i «> t 


7T 7n U 2 + W? g 


zcdpt _|_ i{tjj p —2ujig 


(e ,Wpt + e 


) t) 


’o 

x Tr{dig • G( 

2 f 00 du u 2 e~ uT 


iu ) ■ d * 9 } 


u 2 + w? fl 


J(Up-Ulig)t 


x Tr{djg • [G(iu) - G(-ra)] • d; g } + [ui p -> -ui p \ 

y^r ~^g i f°° d uu 2 (u 2 - UjgUei) 

4^l2ire 0 c 2 J 0 (u 2 + uf )(u 2 + u 2 ei ) 


x Tr{d,;g • G (iu) ■ d~ ig }(e luJat + e 


ioj s t _i_ i(uj s -\-oj e i—uJig)t 


p i(uj s -u) ig )t roo 


du u 2 e uT 


2 J Q (iu — OJig)[iu — LOei) 

x Trjdie • [G(rn) - G(-m)] • d^} 

e i(u,+ Uei )t roo du u 2 e -uT 


+ ~ 


TT J Q U 2 + UJ 2 t 
x Trjdie • G(iu) • d- ie } 


^ e iui s t _|_ e i(w s —2w ie )t\ 


2 I" 00 du u 2 e~ uT 


Jo u 2 + w? e 


,i(ui s —Uie)t 


x Trjdje • [G(m) - G(-iu)] • d- e } + [w s ->• -w s ] 

~^ P gi f°° d UU 2 (U 2 - UigLUei) 


n 


27re 0 c 2 Jq ( u 2 + w 2 )(u 2 + w 2 J 


x Tr{djg • G(iu) • d~ ie }(e iuipt + 


du u 2 e uT 


Jo (iu + UJig )(iu + to ei) 
x Trjdjg • [G(iu) - G (-iu)} ■ d^J 


x [e 


i(uJp—UJig)t I p i{Wp+OJei)t 


]/2 


7 r 


Tr{djg • [<G(w ig )e 1 11 ^ 


i(UJp+CJ ei — UJig)t 


w xg wet 

2 (r*( 

J ig 

,2 i 




i(tJp+Wei— U)ig)t 


0 (iu + uJig)(iu + u ei ) 


x Tr{d ie - • [G(itt) - G(-iu)] • d- g } + [u s -»■ -w s ] sta tes, = w* Vl 


— w eiG(w e j)e*“ p ] ■ d^J +[w p —»• —w p ]. 

As in Sec. IIII Cl the tilded states |i), |g) and \e) belong to 
the same energy levels as the states i, g and e respectively. 
Equal energies have been assumed for all the intermediate 
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